Preprint typeset in JHEP style - HYPER VERSION 



arXiv: 1004.2055 [hep-th] 
DAMTP-2010-28 



Holographic studies of quasi-topological gravity 



Robert C. Myers," Miguel F. Paulos^ and Aninda Sinha" 

° Perimeter Institute for Theoretical Physics, Waterloo, Ontario N2L 2Y5, Canada 

^ Department of Applied Mathematics and Theoretical Physics, Cambridge CBS OWA, UK 



E-mail: rmyers, asinha@perimeterinstitute.ca, m.f.paulos@dciintp.cam.ac.uk 



Abstract: Quasi-topological gravity is a new gravitational theory including curvature-cubed 
interactions and for which exact black hole solutions were constructed. In a holographic 
framework, classical quasi-topological gravity can be thought to be dual to the large Nc limit 
of some non-supersymmetric but conformal gauge theory. We establish various elements of 
the AdS/CFT dictionary for this duality. This allows us to infer physical constraints on the 
couplings in the gravitational theory. Further we use holography to investigate hydrodynamic 
aspects of the dual gauge theory. In particular, we find that the minimum value of the shear- 
viscosity-to-entropy-density ratio for this model is r//s ~ 0.4140/(47r). 
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1. Introduction 

The AdS/CFT correspondence has proven a fertile ground for investigating the properties of 
strongly coupled gauge theories 0, ^, in particular the thermodynamic and hydrodynamic 
properties of these gauge theories at finite temperature |^. However, such investigations 
face acute limitations because at present, we have an insufficient understanding of string 
theory in interesting holographic backgrounds, i.e., in spacetimes with Ramond-Ramond 
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fields. Hence the examination of holographic gauge theories is primarily confined to both 
the limit of a large 't Hooft coupling A and a large number of colours N^. where the dual 
gravitational theory corresponds to (semi-)classical Einstein gravity with a two-derivative 
bulk action. However, it is understood that accounting for higher curvature interactions, or 
higher derivative interactions more generally, within a perturbative framework allows one to 
begin to consider finite A and finite corrections ||5|, ^ ^, ^, ^, An alternative point 
of view would be that admitting such higher curvature (or higher derivative) interactions 
introduces new couplings amongst the operators in the dual CFT, thereby broadening the 



universality class of dual CFT's which one can study with holography ||llj, |12|, y^. If one 
examines a point in the space of CFT's where these new couplings are finite, the higher 
curvature terms will now make finite contributions in the analysis of the dual gravity theory. 
However, if any higher curvature term were to become important, the normal expectation 
is that an infinite number of such terms will become important at the same time as the 
background curvature must have reached the string or Planck scales. The relevance of all 
these terms is really signalling that one has entered a regime where the dual gravitational 
theory cannot be described by a local quantum field theory. Hence to properly investigate the 
effects of these finite CFT couplings, one is brought back to the question of understanding 
string theory in interesting holographic backgrounds. 

However, a traditional avenue to progress in theoretical physics is the study of simplified 
or toy models which might provide insight into the behaviour of some complex physical 
system of interest. Recent work with Gauss-Bonnet (GB) gravity showed that the utility of 
such toy models in a holographic framework |13, 14, [l^, 17, 15]. In this case, the usual 
Einstein action is supplemented by a certain curvature-squared interaction, which corresponds 
precisely to the four-dimensional Euler density. With this extension of the usual Einstein 
action in the five-dimensional bulk gravity theory, the class of holographic models is extended 
to allow independent values of the two central charges a and c of the dual CFT 2C]. 
Further it was found that GB gravity still captures certain fundamental constraints which 
can also be inferred from direct considerations of CFT's alone. In particular, consistency of 
the CFT constrains the central charges to obey fTl]] : 1/2 < a/c < 3/2. Hence GB gravity (or 
more generally Lovelock gravity in higher dimensions |2^]) provides an interesting toy model 
to examine questions related to holographic hydrodynamics, or perhaps the holographic c- 
theorem |22, 

Motivated by the success of holographic studies of GB gravity, this holographic model 
was recently extended with the introduction of a new curvature-cubed interaction in quasi- 
topological gravity [^]. The progress with GB gravity relies on the fact that even though 
this is a higher curvature theory of gravity, the holographic calculations in this model are still 
under control. This control stems from two properties of the theory: the equations of motion 
are only second order in derivatives and exact black hole solutions have been constructed. 
In quasi-topological gravity, exact black hole solutions are again readily constructed but 
on general backgrounds the equations of motion will be fourth order in derivatives |25]. 
Remarkably, however, the linearized equations of motion in an AdSs background are again 
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second order and in fact, match precisely the hnearized equations of the Einstein theory 
[p^ ]. As we will show in the following, these properties are sufficient to allow us to examine 
many interesting features of the holographic framework established by this new toy model. 
The new curvature-cubed interactions again expand the class of CFT's which can be realized 
with this model. In particular, the new couplings are generalized such that the dual CFT 
will not be super symmetric |11] and so this holographic model may provide new insights on 
non-supersymmetric gauge theories with a conformal fixed point. 

One aspect which we examine with this new holographic model are the hydrodynamic 
transport coefficients of the dual CFT, in particular the shear viscosity. It has been observed 
that the ratio of shear viscosity to density entropy of typical holographic fluids is extremely 
small in comparison to typical fluids for ordinary matter . Originally it was conjectured that 
these holographic calculations provided a universal lower bound, namely rj/s > I/Att. One 
piece of circumstantial evidence for this conjecture came from string theory were it was found 
that the leading finite A corrections always raised r//s above the bound for supersymmetric 
plasmas where c = a IH, ^, However, it is now understood this KSS bound can be 
violated in string theory duals of plasmas where c 7^ a by the effect of new curvature-squared 
interactions in the gravitational action [^, |l^. However, the string theory constructions where 
these higher curvature terms are under control only allow for small perturbative violations 
of the KSS bound. General arguments still suggest that the ratio of the shear viscosity to 
entropy density should satisfy some lower bound ^] and so the question naturally arises 
as to the precise nature of such a bound. Hence it is certainly of interest to explore situations 
where finite violations of the bound occur and the toy models above provide a framework 
for such explorations. In particular, it would be interesting if one was able to show that r//s 
could be pushed to zero without producing any other pathologies developing in the theory. 
With quasi-top ological gravity, we find that r]/ s reaches a non-zero lower value in a particular 
corner of the allowed space of gravitational couplings. 

An outline of the rest of the paper is as follows: In section ^, we present a brief review 
of quasi-topological gravity in five dimensions and the black hole solutions of the theory. We 
begin to establish the AdS/CFT dictionary for this gravitational theory in section ^ with 
a calculation of the central charges of the dual CFT. In section ^, we adapt the scattering 
experiments in the CFT of |11] to a holographic calculation. These computations yield directly 
the fiux coefficients t2 and t4, but combined with the expressions for the central charges, we 
are also able to express the coefficients B and C, which determine the three-point functions 
of the stress tensor, in terms of the gravitational couplings. Next in section ^, we consider 
various constraints on the gravitational couplings which are required to ensure the physical 
consistency of the dual CFT. We consider three independent constraints: positivity of the 
central charge c, positivity of the energy fiuxes in section Q and avoiding violations of causality. 
In section ^, we examine the hydrodynamic behaviour of the CFT plasma. In particular, we 
find that the minimum value of the ratio of the shear viscosity to the entropy density in this 
model is r]/ s ~ 0.4140/ (47r). We provide a preliminary analysis of possible instabilities of the 
black holes, or alternatively of a uniform plasma in the dual CFT at finite temperature in 
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section |^ We conclude with a brief discussion of our results and future directions in section 



2. Review of quasi-topological gravity 

We begin with a review of some salient features of quasi-topological gravity. We focus on 
the five-dimensional version of the gravity theory, which would be dual to a four-dimensional 
CFT. The bulk gravity action can be written as p5[|: 



1 

2a 



where <^4 is the four-dimensional Euler density, as used in GB gravity 

^4 — R^vpaR^ ^ ^Rpu 



(2.1) 



Xa — RavoaR^'^'^^ — ^RauR^'^ + R? i (2.2) 



and Z'^ is the new curvature-cubed interaction 

Z'^ = + ^ (21 R^.p^R'^^P'^ R - 120 R^.p^R^^" ^R''^ (2.3) 

+UAR^,yp„R^'PR'''' + 128 R^^R^PRp^ - 108 R^" R^^ R + 11 R^) . 

The AdS vacua of this theory have a curvature scale given by 



L2 



(2.4) 



where the constant foo is determined as one of the roots of 

1 - /oo + A/i + m/^ = . (2.5) 
Note for any choice of the couplings A and /i, there is at most one ghost-free AdS vacuum 



which supports nonsingular black hole solutions, as described in detail in [25] . The solutions 
describing planar AdS black holes take the form 

ds^ = 4 f -4^ dt^ + dxl + dxl + dxi] + dr2 , (2.6) 



where /(r) is determined by roots of the following cubic equation: 

l-/(r-) + A/(r)2 + ^/(r)3 = 4 . (2.7) 

For the relevant solutions, the black hole horizon occurs at r = tq, which is easily seen to 
yield /(r = ro) = as a solution of the above equation. The Hawking temperature is given 
by 



7rL2 /.i/s ■ 

Joo 



(2.8) 
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The energy and entropy densities are simply calculated as |2 

Further note that these relations satisfy p = jTs, as expected for a four-dimensional CFT 
(in the absence of a chemical potential). 

Apart from finding exact black hole solutions, another remarkable property of quasi- 
topological gravity is that the linearized graviton equations in the five-dimensional AdS vac- 
uum take the form |p5|| : 

V^/lab + VaVfe he' - VaV'Kb - ^b^'h^a (2.10) 



1 '1 - 2A/oo - 3m/^) 



2 



L2 L2 

Here (7^"^ is the background AdSs metric and L is the curvature scale in eq. ( |2.4j ). Hence the 
observation is that in the AdSs background, the linearized equations are only second order 



in derivatives. In fact, up to an overall factor, the above equations (2.10) are precisely the 



same as the linearized equations for Einstein gravity in an AdSs background — for example. 



see 1 27, 28 1. This result contrasts with that for a generic B?" action which would yield fourth 
order equations of motion.^ However, the same occurs here for quasi-top ological gravity on 
a general spacetime geometry. That is, on general backgrounds, the linearized equations are 
fourth order in derivatives for the present theory as well. 

3. Central charges 

In this section and the following section, we develop the dictionary relating the couplings 
in five-dimensional quasi-topological gravity to parameters which characterize the dual four- 
dimensional CFT. Since we are only dealing with the gravitational sector of the AdS theory, 
we are looking to examine the behaviour of the stress energy tensor of the CFT. Two such 
parameters are the central charges, c and a, of the CFT. We calculate these through their 



appearance in the trace anomaly |1T^, using the now standard holographic approach of |31]. 



The central charge c also fixes the coefficient of the leading singularity in the operator product 



of the stress tensor with itself [P3, 33, p4]. Hence as a verification of our first calculation, we 



also determine c from examining the two-point function in section |3.2| . 
3.1 Holographic trace anomaly 

The two central charges of a four-dimensional CFT can be defined by the trace anomaly that 



arises when the CFT is placed on a curved background geometry |19|: 



(T^'^) = " ;t'4, (3.1) 

^ ' 167r2 167r2 ^' ^ ' 



^Recently, other theories of curvature-cubed gravity with exceptional properties were identified in |3o|] . 
In fact, up to a contribution proportional to the six-dimensional Euler density, the curvature-cubed interaction 
constructed in five dimensions by [B0| is identical to that studied here. 
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where X4 is the four-dimensional Euler density, whose structure is given in eq. (|2.2D (although 
here, is evaluated for the four-dimensional background metric of the CFT), and I4 is the 
square of the Weyl tensor, i.e., 



J4 — '^abcd ^ — n,abcd ^ — ^ ^ab^ 3 ' 



(3.2) 



In order to compute c and a for the CFT dual to quasi-topological gravity, we follow the 
holographic procedure described in [^]. We should note that modifications to the central 



charges from interactions were examined previously in |2C] while perturbative corrections 
coming from interactions were considered in [^^. Efficient methods, i.e., 'short cuts' to 
calculate the holographic trace anomalies for an arbitrary gravitational action are discussed 



in H, |4|. 



Following 1 31 1, we begin with the Feffer man- Graham expansion 



J^d/J2 + Mdx«dx^ 



(3.3) 



with 



9ab = g{0)ab + P9{l)ab + P 9{2)ab H , (3.4) 

where the boundary metric g(^o)ab corresponds to the background geometry of the dual CFT. 



The next step is to substitute this expansion of the metric into the gravity action (2.1). On- 
shell 5((2) drops out and we are left with an action involving ^r^g) and g^^iy To extract the 
conformal anomaly, we focus on terms which when integrated produce a log divergence. This 
leads to 



/5(0) 



abed 



(3.5) 



-A R^'^'^'g^Dab + B tr 5(1) + C tr gf,^ + D (tr 5(1) ) 



where, e.g., R^^^ corresponds to the Ricci tensor calculated for the boundary metric 5(0)06 • 
Further e defines a UV regulator surface which cuts off the radial integral. The constant 
coefficients appearing in the above expression (^^) are given by 



h 

A 
D 



3/^/i) 



-(I-2A/00 
1 



6 



(7 - A/00 + 5pf^ 



At2 = ts , 
) = -2B 



L3 



(3.6) 



/73 f 3/2 ■ 



Next we eliminate g(i)ij using its equation of motion and then following |31], we can interpret 
the result as 

/:^-loge- j d'^x^iTa^). (3.7) 
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Hence comparing the coefficients of the various terms involving the background curvatures 
with eq. (|3.lD, we find 



^3 I 

c = vr^-Tj-T^ (l - 2A/oo - 3///^) , (3i 



li 1 

a = ^^-^^{:^-^^foo + ^^lflo) . (3.9) 



"P /oo 

Given these results in eqs. ( |3.8| ) and (p^), we also have 

c-a 4/oo(A - 3/i/oo) 



c 1 - 2A/oo - 3/i/, 

3.2 Two-point function 



(3.10) 



Now we turn to computing the two-point function of the stress tensor as an alternative ap- 
proach to determining the central charge c. It is known |3^, 33, 36] that in a four-dimensional 
CFT2 ^ 

{Tab{x)Tcd{x')) = j^^——^lab,cd{^-^') (3-11) 

where 

1 1 

'^ab,cd{^) = 2 (^<ic{x)Im{x) + Iad{x)hc{x)) - -r]abVcd 

and Iab{x) = rjab - 2^^^ . (3.12) 

This structure is completely dictated by the constraints imposed by conformal symmetry 
[ p^ , |36[ . The coefficient Ct is related to the central charge c which appears as the coefficient 
of the (Weyl)^ term in the trace anomaly ( p.l[ ): 

40 

Ct = —c. (3.13) 

In order to compute Ct, it is sufficient to focus on the specific case {T^yTxy). To 
determine this two-point function, we will turn on a perturbation r'^hxy{r, z)/L?' in the AdSs 



background, i.e., in eq. ( |2.6| ) after setting rg = 0. The quadratic action for h^y = (I) can be 
written as 



where 



/2 = ^ y d^'x {Kr{dr(t>? + K,{d,(t>f + drT) , (3.14) 



».5 .1/2 

i^r = (1 - 2A/oo - 3/z/^) , (3.15) 

= YK- (1 - 2V°o - 3/i/i) . (3.16) 



We assume a Minkowski signature for the metric. Hence in eq. (3.12), Xa ~ rjabX (i.e., xq = —t). 
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The details of T are unimportant since the drT contribution is canceled by a generalized 
Gibbons-Hawking boundary term S. Upon making the ansatz 



the equation of motion for (p reduces to 



The general solution can be written as 



(3.17) 



(3.18) 



(3.19) 



where l2{x) and K2{x) are modified Bessel functions of the first and second kind, respectively. 
In order to fix Hy(r = oo) = 1, we set Ci = p'^L^/ (2foo) and C2 = 0. Using the equations of 
motion, eq. (3.14) can then be rewritten as 



/2 = d^xdr{Kr(t>dr 



(3.20) 



Substituting in the solution, we make a Fourier transform and extract the term proportional 
to log \p\ , which amounts to ignoring all contact terms: 



( T^y Txy ) (p) 



L3 



3/2 
pjoo 



:i -2A/00- 3^/4)/ log 



(3.21) 



In order to compare the above result with eq. ( p.ll ), it convenient to recast the latter in 

(3.22) 



Fourier space following [37|. We re-express the two-point function as |36| 

1 



{Tab{x) Tcdix' 

where the tensor S'-' satisfies 



40 



{x - x')^ ' 



^abef ,cdgh k'^k^k^k^ = 777 ( 2 kak^kck^ — —k^ikakcTjlfd + k^kfTjad + kakdTJbc + ^b^d^ac) (3.23) 



24 



+k {kakf,ricd + kckdVab) + 7;ik ) {Vacribd + VadVbc) - {k ) VabVcd ■ 



In the case of interest with (ab) = {xy) = (cd), the factor involving S*-" and the four derivatives 
simply evaluates to 

//16. (3.24) 
Hence in momentum space, the two-point function can be written as |37] 



{T^yT^y){p) = ^P'' i d^x 



P log IpI + (analytic in p) . (3.25) 
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Comparing eqs. ( 3.21 ) and ( ^.25 ), we find 

40 1 



Ct = ^ 



7r2 f3/2 



(1 - 2A/oo - Sfifl) 



(3.26) 



and finally using eq. ( 3.13|) , we have 

,L3 1 



C = TT 



a .3/2 



^ (1 - 2A/oo - 3^/^) 

P foo 



(3.27) 



This expression precisely matches that in eq. ( |3.8D which was found using the holographic 
trace anomaly in the previous section. 



4. Holographic computation of energy fluxes 



At this point, our holographic dictionary contains two entries. That is, eqs. ( p.8| ) and (|3.9| ) re- 
lating the central charges of the four-dimensional CFT to the couplings of the five-dimensional 
bulk gravity theory. However, the quasi-topological gravity ( |2.lD is characterized by three in- 
dependent dimensionless parameters: A, /x and L/ip. Hence we need to extend the dictionary 
further by identifying additional parameters which play an analogous universal role in the 
dual CFT. Further, since we are only dealing with the gravitational sector of the AdSs theory, 
we must find parameters governing the behaviour of the stress tensor in the CFT. A natural 
next step is to consider the three-point function of Tij, as was extensively studied by [33, 36]. 



There it was shown that conformal symmetry and energy conservation are powerful enough 
to determine the the three-point function up to three constants, which are labeled A, B and 
C in [^]. In fact, the two central charges can be expressed in terms of these three param- 
eters, as we will elucidate below — see eqs. ( 4.35 ) and ( 4.36| ). Further as discussed in |11], 
constructing a holographic model which can explore the full range of these CFT parameters 
requires the introduction of curvature-squared and curvature-cubed interactions in the bulk 
gravity theory. In fact, this was the primary motivation for constructing the quasi-topological 
gravity theory. 

One can perform a holographic calculation of the full three-point function [^], however, 
extending these calculations to quasi-topological gravity proves to be extremely challenging. 
Therefore, we choose an indirect route to determining these coefficients here. In particular, 
we construct a holographic description of a particular thought experiment proposed for four- 
dimensional CFT's in The experiment consists of first producing a disturbance, which 
is localized and injects a fixed energy, with an insertion of the stress tensor £ij T^^ , where £ij 
is a constant polarization tensor. Then one measures the energy fiux escaping to null infinity 
in the direction indicated by the unit vector n. The final result takes the form 



E 
An 



l + t2 



+ t4 



EijU n- 



i|2 



(4.1) 
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where E is the total energy. The structure of this expression is completely dictated by the 
symmetry of the construction. Hence the two constant coefficients, t2 and t^, are parameters 
that characterize the underlying CFT. Our holographic computation of t2 and ^4 will extend 
our AdS/CFT dictionary to the point where the three independent gravitational couplings 
will be related to three independent parameters in the dual CFT. 

Note that the (negative) constants appearing in eq. ( [4 .11 ) in the two factors multiplied 
by t2 and were chosen so that these factors contribute zero net flux when integrated over 
all directions. The negative sign of these constants leads to interesting constraints on the 



coefficients t2 and 14, which we discuss below in section ^ 
4.1 Field theory calculations 



Let us first consider the discussion of |11| in more detail. Again, we begin by making a small 



localized perturbation of the CFT in Minkowski space with metric ds'^ = —dt^ + Sijdx'^dxK 
With time this perturbation spreads out, and sufficiently far away one may imagine taking 
successively larger concentric two-spheres through which one is measuring the energy flux. 
Let us parameterize the points on these two spheres by a radius r (i.e., = XjX* as usual) 
and a unit vector n = {n} ,v? ,n^). Then the energy flux measured in the direction given by 
n is given by 

/+00 
dt TK,{t,rn) n\ (4.2) 
-00 

As it stands, the flux in eq. ( |4.2| ) includes contributions from both the past and future 
null boundaries of Minkowski space. To separate out the future contribution which we are 
interested in, we select one of the coordinates and construct light-cone coordinates x^ = 
tibx^. Then it is clear that the integral above has two contributions, namely one from future 
null infinity x^ — t- +00 and another one from past null infinity x~ — t- —00. We will take only 
the former.'^ For large r, it is convenient to write 

^2^(^+_^-)2^(^l)2^(^2)2 -+^00 (^+)2^i^(yl)2^(y2)2)^ ^^^^^ ^1,2 ^ ^l,2/^+_ 

Therefore we obtain 



,n 



(4.3) 



/ + 00 
dx' [r+,(x+,x-,n) + r_,(x+, 
-oo 

Motivated by the preceding, we define new coordinates y": 

+ 1 - - (Xl)2 + (x2)2 ^2 

y^ = — T, y =x , , y 4.4 

In terms of these coordinates, the desired energy flux is measured = 0. Further it is not 
difficult to show that on this surface, we have 

y''' = • (4-5) 

1 + n'^ 

^Note that the following discussion overlooks the flux contribution on the 'hemisphere' at x~ +oo. 
However, our primary concern is the functional dependence of S{n) and so this does not affect our results. 
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Now we transform the energy momentum tensor from x"" to y'^ coordinates, as usual 

^dy^dy^ y 
dx^ dx^ ^'^ ■ 



(4.6) 



Here we are simplifying our notation with the superscripts, x and y, to indicate in which 



coordinate system the stress tensor is written, i.e., T^__ 



-y y 



Now at y"*" = 0, we obtain 



++ 



(4.7) 



Tl_ = {{y^f + (y2)2) Tl^ , T^l_ = 



- 1,2 



2yi,2 Ty__. 



Hence we see that on this surface {i.e., y^ = 0), there is a single relevant component of the 
energy momentum tensor in y coordinates. Using these relations, we rewrite eq. (|4.3|) as: 



£{n) = 



, _ Ty^^iy^ = 0,y-,y\y^) 

-oo m-^ 



(4. 



with n = 2/(1 + n^). 

Further, we note that transforming from the x° to y"" coordinates produces a Weyl scal- 
ing of the metric. Hence it is natural to use this transformation to perform the conformal 
transformation: 



-dy+dy- + (dyi)2 + {dy^f 



ds2 = -dx+dx- + (dx^ + (dx'^f ---^ 

-2- -dy+dy- + (dyi)2 + (dy2)2. 



ds' 



Then the energy momentum tensor transforms 



Tab 



dx 



dy 



^^^"^ dx^ dx'^ 

Qya Qyb 



Ted ■ 



(4.9) 



(4.10) 



In particular, we have = r^_/(y+)2 and therefore eq. (4^) becomes: 



f(n) = 



+ 00 



dy Ty__{y^ = 0,y ,y ,y ). 



(4.11) 



Following [11 1, we wish to consider the expectation value of this flux operator £{n) for a 
particular state 

{0\Oi£{n) Oe|0) 



(^(n)) 



(0| Ol, |0) 



(4.12) 



In the present discussion, the operator is assumed to be a localized insertion of the stress 
tensor of the form 

(4.13) 



Oe = j d'^x e'^ Tij e-'^^ ^{x/a) . 



Here iplx/a) is some profile which localizes the insertion to x'^ = on the scale a. We assume 
E ^ 1/a and so the energy of the insertion is E up to order l/cj corrections. Finally since 
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the stress tensor enters this construction, the operator also contains a constant polarization 
tensor eij which we assume only has spatial components. The symmetry of this construction 
then dictates that the flux expectation value ( [4. 121) takes the form given in eq. (4J). Further, 
it is clear that by construction the result is completely determined by the three-point function 
of the stress tensor. Hence the flux parameters t2 and ^4 appearing in eq. ( [4.1| ) will be related 
to A, B and C, the parameters controlling this three-point function. 

4.2 Holographic calculations 

The x"" and coordinates defined in the previous section are easily extended into the AdSs 
bulk with 

ds2 = _ (-dx+dx- + idx^f + (dx^f + dz^) (4.14) 

= ^ (-dy+dy- + (dyi)2 + {dy^f + du") . (4.15) 

Recall that L is the curvature of the AdSs geometry, as defined in eq. ( p.4p . To relate these 
two coordinate systems, it is convenient to describe AdSs as the hyperbola 

in a six-dimensional Minkowski space with (-,-,+,+,+,+) signature. Note that we reach the 
boundary of AdSs by taking W^^ large. The previous coordinates are mapped to the 
coordinates with 

W-^ + W'^ = —, W'' = L— for = 0,1,2,3, (4.17) 

z z 

W^ + W^ = -, W-^ + W^ = -P — , W-^-W^ = -—, W^'^ = L — . 

u u u u 

Note that z and u are mapped to two orthogonal null surfaces in the W^'^ space. Further the 
powers of L are slightly different in the second line above to ensure that the (engineering) 
dimension of the coordinates is properly accounted for, i.e., u is dimensionless while has 
dimensions length~^. With eq.( 4.17 ), we can relate the {x°',z) and {y'^,u) coordinate systems 
m eqs. (PI ) and (|1|) as 

I 1 1 9 X ' z 

x+ x+ x+ 

x+ x+ 

Notice that on the asymptotic boundary z = 0, the above coordinate transformation reduces 
to that given in eq. ([4.4|). Further with y+ = and any finite value of u, we are on the AdSs 
horizon at 2; = 00 in the (x", z) coordinates. 

As commented above, in calculating the flux expectation value in eq. ( [4.12 ), we are essen- 
tially determining a specific component of the three-point function of the stress tensor. Hence 
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in our holographic description, we must first introduce appropriate metric perturbations /i^jy 
in the AdSs bulk which couple to the dual insertions of Tab- We then evaluate the on-shell 
contribution to the cubic effective action for the graviton insertions. 

As discussed in |25], in general, the equations of motion for quasi-topological gravity 
involve higher derivatives. Hence one would expect that linearized equations of motion for 
the metric perturbations here are also higher order. However, it was observed in |25| that in 
fact these linearized equations for gravitons propagating in the AdSs vacuum match precisely 
the second order equations of Einstein's theory, up to some overall (constant) coefficient, as 
shown in eq. ( |2.1Cl| ). This makes the following calculations much simpler as we may borrow 
previous results |lT| , |2^ for the graviton solutions in Einstein gravity. Hence while the higher 
derivative contributions in quasi-topological gravity are essential to producing a nonvanishing 
value for t4 in the dual CFT, they only contribute through the three-point interactions in the 
following. 

We first consider the flux operator £{n) in eq. (4.11). It is natural to use the (y", u) coor- 
dinates in eq. ( [1.15D and the standard AdS/CFT dictionary advises us that (y) couples 

to hj^^{y'^,u = 0). Considering first a localized insertion /i+_|_(y",n = 0) = 6{y^)6{y'^)d{y~^), 
the bulk solution is given by 

h,,iy-, y-,y\y\ u) = + (,i)2 + (,2)2)4 • (4-19) 

As noted above, we are using the same solution here as in [|ll| ] because the linearized equations 
of motion for perturbations around AdS^ in quasi-topological gravity are the same for Einstein 
gravity [^]. To obtain the operator if(n), we then integrate in y~, as well as multiplying by 
an overall factor of 17'^ and performing a translation in and y^^ to obtain 



h++{y^,y^ 



8d{y^ 



u 



(1 + n^f (n2 + (yl - y'^f + (y2 _ y/2)2)3 ' 



(4.20) 



with y' 



'1,2 



n 



1,2 



/(I + n'^), as in eq. (| 



In fact we can make this insertion at a nonlinear level in the bulk gravity theory, following 
n and |l|]. To achieve this, we consider the shockwave background: 

Z2 



ds' = -J 



Ky^) W{y\y^u){dy+f - dy+dy~ + {dy^f + {dy^ f + du' 



(4.21) 



This metric solves the full equations of motion coming from eq. ( |2.1| ) provided that yV{y^,y'^,u) 
satisfies the equation of motion 



dlW-- duW + d^W + d\W = 



u 



(4.22) 



This simple linear equation appears as the equation of motion in Einstein gravity and one 
can readily show that it is not corrected by the higher curvature terms in eq. (2.1) with the 
arguments of 



From our expression for 
W{y\y\u) = %- 



L2 {v? 



in eq. (4.2C), the relevant wavefunction is 

,4 

(4.23) 



u 



+ (yl _y'l)2 + (y2 



72^2^3 
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with y'-*^'^ = + n'^), as before. 

Next we turn to the graviton perturbations dual to the operator insertion Oe in eq. ( [4.13 ). 



To simphfy the discussion, we choose a particular polarization with e^ix^ = 1 = and all 
other components vanishing. Using the {x'^,z) coordinate system in eq. ( 4.14| ), the desired 



operator ([4.13| ) is sourced by a metric perturbation with the boundary value: h^ix2{x°',z — )• 
0) = 2;~^e~*'^^^^'''^ The bulk solution that corresponds to this boundary perturbation is 
then 

= /<^.'.-<f(^-«'-. (424) 

Since the perturbation is completely localised at y'^ = 0, for later purposes, we will 
primarily be interested in the behaviour of hj.ix2{x"', z) on that surface. Following [|ll|, it is 
possible to perform the above integral using the parameterization of AdSs in eq. ( [4.16 ) to 
produce 

KiAW^ 0, W-,W^) e-''''^-/'6Hw') , (4.25) 

where = ± W^. Imphcitly, has been replaced with {W^)^ = 1 — (VK*)^ which 
is the reduction of eq. ( [4.16| ) with = 0. Using eq. ( 4.171) , we may express the coordinate 
dependence in terms of {y°',u) 

K^AV^ ^ 0,y-,y\y\u) ^^J^e^J^y'^^ 6iy') 5{y^) 5{u - I) . (4.26) 



Finally with the coordinate transformation ( [4.18 ), we also transform the tensor indices to 



find that at y^ — t- 0, our metric perturbation becomes 

K-y^ ^ ^ e*""'"^' <^(y') '^(y') - 1) (4-27) 

along with other hy+yi, hy+y2 and hy+y+ components. However, the form of the latter will 
not be important, as we now discuss. Note that the original expression (4.24) was transverse 



and traceless in the {x°',z) coordinates but as a result, the expression produced by simply 
making a coordinate transformation to the (y", u) coordinates is not. However, it is convenient 
to work in this gauge since a great simplification results in the equation of motion for the 
graviton propagating in the AdSs background, i.e., away from the shockwave deformation 
in eq. ( 4.21[ ). Hence at this point, we choose add to eq. ( 4.27[ ) the components required to 



impose transverse traceless gauge in the {y^,u) coordinates. 

The mode above was traceless by construction and so we only need to ensure that the 
transverse condition is satisfied as well, i.e., V^/i^jy = 0. In the present case, the latter can 
be satisfied as long as hyiy2 is accompanied by modes satisfying: 

^y~ ^y^y^ "2 ^y^^y^y^ ' ^y^^y^y^ '2^y^^y^y^ ' 

dy-hy+y+ = ^ [B yl k yl y+ + ^J, 2 /l 2 + ) . (4.28) 
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Together the hyiy2, hy+yi, hy+y2 and hy+y+ components form an independent transverse 
traceless mode. We have verified that with hyiy2(y,u) = LP' /v? (f){y^u) and the remaining 
components chosen to satisfy eq. (|4.2S| ), the equation of motion for (j){y,u) becomes simply 
that of a massless scalar in AdSs (up to interaction terms with the shockwave): 



u 



4:dy+dy- 



0. 



(4.29) 



To find the three-point function, we add these perturbations to the metric ( 4.21 ) and 
evaluate the action ( |2.1| ) on-shell. Then we must extract the terms of the form from 
this result. After integration by parts and using the equations of motion, the cubic effective 
action becomes 



o(3) 



— / d''xV^cpd^_cl)W l-2f^X-3fifl + f^{X-87f^f,)T2 + 2lflfiT^ (4.30) 



where 



U=l, ?/l=j/2=0 



3T2 + 



dtdjW - dudfW - dudjW 



u=l,yl=j/2=o 



(4.31) 



Implicitly here, we are using that with the perturbations given above, i.e., eqs. ( [4.20| ) and 
( 4.26| ), the interaction is entirely localized along y+ = = = and u = 1. Substituting 
the solution ( [4.23| ) for yV{y^,y'^,u) into eq. (4.31), we obtain 

-|- 77-2 



24 



,2 I ^2 



T4 = 180 ( 2 n| - — 



(4.32) 



The expressions above involving and should be interpreted as two independent 50(3)- 
invariant combinations of the unit vector n* and the (implicit) polarization tensor Sij, i.e., 



2 " 4^.. 



2n\n\ 



Sij n n- 



(4.33) 



This was the guiding principle in selecting the two combinations presented in eq. ( [4.31| ). To 
normalize the final result, we must divide by the two-point function {Tyiy2Tyiy2), which is 
essentially the calculation of section 3^. We finally arrive at an expression identical to that 

24/oo(A-87/oo/i) 



in eq. (4T) with 



t2 



1 - 2/ooA - Sflfi 



3780/4/i 
1 - 2/ooA - 3fl^i ■ 



(4.34) 
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4.3 General three-point parameters 

At this point, we return to A, B and C, the parameters in the dual CFT controlling the general 
structure of the three-point function of the stress tensor ||3^, |3^. In fact, these parameters 
fix both the central charges, c and a, and the flux parameters, t2 and t^, in the CFT. Hence 
we can use the results in the previous section and in section ^ to express the three-point 
parameters in terms of the gravitational couplings. 
First, we can express the central charges as [^] 



TT 

480 
7r6 



2880 



{QA-B- IOC) , 
(13^ -2B- 40C) 



Further we have |11, 

15(5^ + m- 12C) 

*2 ~ 



tA 



15(17^ + - 80C) 



(4.35) 
(4.36) 

(4.37) 



9A-B-IQC ' " 4(9^ IOC) 

Given that these four quantities are all determined by the same three parameters, these 
expressions must be redundant. That is, we can see that there is a consistency condition: 

41^ -AB- 20C 



1 4 

-t2 ^ — ^4 

6 45 



6(9^ -B- IOC) 
With the results ( [4.34 ) in the previous section, 

4/oo(A - 3/i/oo) 



1 4 

-t2 H ti 

6 45 



f2 
I oo 



(4.38) 



(4.39) 



1 - 2A/oo - 3/i/2 

Now comparing to eq. ( 3.1C| ) , we see that our holographic results satisfy the required relation 

(H)- 



Combining these expressions (4.35~ 07| ) with the results of the holographic calculations, 
eqs. (4.34), ([4.35| ) and ( [4.36| ), we arrive at the following expressions 

512 L3 1 



A 



B 



C 



97r4 £3 


.3/2 
Joo 


32 L3 


1 


97r4 £3 


.3/2 
Joo 


32 L3 


1 


97r4 £3 


.3/2 
Joo 



foo) 



(1 - 12A/oo + 48^/^ 
(49-318A/oo+4377^/^) 
(23 - I68A/00 + 213^/^) . 



(4.40) 
(4.41) 
(4.42) 



5. Physical constraints 



Having established several interesting entries in the AdS / CFT dictionary for quasi-topological 
gravity, we next consider various constraints on the gravitational couplings that arise to ensure 
the physical consistency of the dual CFT's. We consider three independent constraints in the 
following: 
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5.1 Posit ivity of Ct 

Unitarity of the CFT requires that the central charge Ct or c is positive. This constraint on 
c may seem somewhat mysterious from the point of view of the holographic trace anomaly 



discussed in section 3.1. However, the definition in section |3.2| shows Ct appears in the two- 
point function and so this central charge controls the sign of the norm of CFT states created 
with the stress tensor. Hence given eq. ( p. 26 ), the dual gravity theory must satisfy 



l-2/ooA-3/^M>0. (5.1) 



At this point, we note that multiplying the expression on the left-hand side of eq. (5.1) 



by minus one yields precisely the derivative of the left-hand side of eq. (P^), i.e., the slope of 



the polynomial there evaluated at the root given by eq. (2.5). This comment is related to the 



observation in |25| that this slope appears as a pre-factor in the kinetic term for graviton or in 



the linearized equations of motion ( |2.1C ) in the AdS vacua of the theory. That is, the sign of 



the slope determines whether or not the graviton is a ghost in a particular AdS vacuum — the 
graviton is well-behaved when the slope is negative. Of course, the holographic calculation in 



section 3.2 shows that Ct is precisely determined by the graviton propagator and so it is no 
surprise that the same factor appears in both places. Further, we note that in AdS vacua with 
a ghost-like graviton, this pathology would make a prominent appearance as non-unitarity in 
the dual CFT since Ct would be negative. The analogous observations were made for GB 



gravity in |17 



5.2 Positivity of energy fluxes 

Turning to the expression of the energy flux in eq. ( [4.1D , we note that the two factors multiplied 
by t2 and i4 were normalized to give a vanishing contribution to the net flux when integrated 
over all directions. Hence depending on the specific direction, these factors may give either a 
positive or negative contribution to (<S(n)). Further, it is easy to see that if the coefficients t2 
and t4 become too large, the energy flux measured in various directions will become negative. 
Following [^], avoiding this problem then imposes various constraints on these coefficients:^ 

12 , , 

Tensor : 1 - -t2 - — ^4 > , (5.2) 
6 15 

1 2 

Vector : 1 + -t2 - — ^4 > , (5.3) 

15 

1 8 

Scalar : 1 + -t2 H ^4 > . (5.4) 

3 15 

If the polarization tensor is chosen with e^i^'^. = 1 = Ex'^x^ all other components vanishing, 
as in section |4.2| , the tensor, vector and scalar constraints above correspond to demanding 
a positive flux with n| = 1, = 1 (or n\ = 1), and n\ = 1/2 = n\, respectively. Using 



complementary analysis in ||39|] produced a constraint equivalent to eq. (5.4), again as a positivity 
constraint on the three-point function of the stress tensor. 
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eq. ( 4.34 ), we translate these constraints on t2 and i4 to constraints on the gravitational 
couplings: 



Tensor : 1 - 10/ooA + 189/^/i > , 
Vector : 1 + 2/ooA - 855/^^ > , 
Scalar: 1 + 6/ooA + 1317/^^ > . 



(5.5) 
(5.6) 
(5.7) 



In the present case, these constraints confine the higher curvature couplings of quasi-topological 
gravity to lie within a small region in the (A, ^)-plane, as shown in figure ||. 



^ 0.000 




Figure 1: The allowed region in (A,/i)-plane satisfying the constraints appearing in eqs. dsiS-js^). 
Within this region, the energy flux ( [f.l[ ) in the dual CFT is positive for any direction. The (blue) 
segment on the A-axis within the allowed region matches precisely the allowed values of the coupling 
in five-dimensional GB gravity. 



Setting fi 
straints from eqs. (15. 



reduces the theory to GB gravity and one recovers the expected con- 
]^) in this limit [[ 



12|]. First with /x = 0, eq. (|2j) yields f, 



^ (l — \/l — 4A) for the ghost-free AdS vacuum. Then, for example, the tensor constraint 



( |5.5D reduces to 5\/l — 4A — 4 > or A < 9/100. Similarly, the vector ( |5.6| ) and scalar ( |5.7| ) 
constraints yield A > —3/4 and A > —7/36, respectively. Hence to maintain positive energy 
fluxes in all directions, the curvature-squared coupling must lie in the range — ^ < A < as 
expected |jl5|, While the latter combines the results for the tensor and scalar constraints, 
the inequality arising from the vector constraint also matches the previously derived result 
|l2| . The allowed GB theories are illustrated in figure || as the blue segment on the A-axis 
(i.e., 1-1 = 0) within the allowed region. 
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5.3 Causality constraints 

Constraining the gravitational couplings by demanding that the dual CFT respects causality 
was first explored in the context of five-dimensional Gauss-Bonnet gravity [0, 13 1. In this 
analysis, one considers graviton fluctuations that probe the bulk geometry. The AdSs vacuum 
{i.e., eq. (2.6) with /(r) = foo) is Lorentz invariant in the CFT directions and so no violations 
of causality would be found with this bulk spacetime. Instead, the black hole solution provides 
a background where Lorentz invariance is broken and in certain instances, the dual CFT 
plasma supports superluminal signals. Hence one constrains the gravitational couplings to 
avoid the appearance of such superluminal modes. The original analysis |14, ^] of GB 
gravity only considered gravitons polarized transversely to the momentum direction, in what 
is conventionally called the tensor channel. The analysis was later extended to the shear 



and sound channels in [15, O. These causality constraints have since been extended to GB 



gravity in higher dimensions [17, 16, 18| and more generally to higher order Lovelock theories 

In all of these cases, it was found that the causality constraints precisely match those 
arising from requiring positive energy fluxes. In particular, for five-dimensional GB theory. 



these constraints are exactly equivalent to those presented in eqs. (|5.5H5.7D with /i = 0. 
However, it has been shown that this matching does not appear in general, in particular 
for cases where the gravitational equations of motion are not second order [^]. Hence, in 
general, one has two independent sets of constraints, one required by the positive fluxes and 
a second determined by the absence of superluminal modes. In quasi-topological gravity, the 
linearized equations in a general background (and as we will see, in a black hole background) 
are fourth order in derivatives and so we do not expect that the previous constraints ( ^.5| - |5.7D 
will be reproduced by the causality analysis. However, unfortunately our final results here 
will be similar to those in [12[. That is, we find no evidence of causality violation once the 
curvature-cubed coupling fi is turned on in quasi-topological gravity. 

There is a broad literature discussing causality in general field theories [40[. The key 
property characterizing how quickly signals propagate is the speed with which a wave-front 
propagates out from a discontinuity in some initial data. This front velocity is given by 

vf'''^'= lim Re{u)/q. (5.8) 

|g|-i-oo 

That is, we are interested in the phase velocity of modes in the limit of infinitesimally short 
wavelengths. Hence in a relativistic field theory, we would require that v^'^°^^ < 1 in order to 
avoid any acausal behaviour. 

In the present holographic framework, we need to determine the front velocity of signals in 
the dual CFT. That is, we must determine this velocity for excitations dual to various graviton 
channels in the bulk spacetime. Consider the black hole background given in eq. ( |2.6| ) and 
define the new coordinate p = r^jr^ ■ The metric becomes 

= ^ (-M + dx? + dx2 + Axl] + dp2 . (5.9) 



- 19 - 



For simplicity, we will focus on tensor perturbations of the form 

„2 



= -^e--*+^'?-=<A(p) (5.10) 



propagating on the background given in eq. ( |5.9D . The determination of the front velocity 
for such modes was described in detail in |jl3|. For our present purposes, it suffices to derive 
the effective speed in the CFT directions at large momentum and frequency by focusing on 
the contributions coming from the t and z derivatives in the linearized equations of motion. 
That is, the full linearized equation for (pip) takes the form 



d, 



The function C^^^ contains two terms, one independent of q and the other proportional to g^. 
In contrast, the C^^^ function is a sum of terms proportional to uP' ^q^ ^uP'q^ and q^. In the 
large momentum and frequency limit, the radial derivatives can be neglected and essentially 
only the term is relevant above. By setting this term to zero, one finds an effective 
'dispersion relation' relating the frequency to the momentum. As the final result is quite a 
complicated expression, let us approach it in several steps. 

We start by setting both A and /i to zero, in which case our theory reduces to Einstein 
gravity. Then the dispersion relation is simply 

= c^2_Mg2_ (5.12) 

/oo 

For Einstein gravity, /(p) = \ — (and /oo = !)• Therefore the pre-factor multiplying 
q^ above is less than one for any finite radius and we expect that these excitations always 
propagate at less than the speed of light, i.e., /q^ < 1. Next we consider GB gravity with 
A ^ and /i = 0. In this case the dispersion relation becomes 

= a;2 (1 - 2A/(p) + IpXfip)) - M q^l - 2\f{p) + IXpfip) - Ap^Xfip)) (5.13) 

Joo 

Using the GB black hole solutions and expanding the above expression near the AdS boundary 
yields 



Hence we expect that preventing a superluminal front velocity requires 1 — 10/ooA > 0, which 
matches precisely the tensor constraint ( ^.5| ) with p = 0. Of course, this agreement for GB 
gravity was previously noted |11, 15, O. 



Now for the full theory with both A 7^ and p ^ 0, the story is quite different. As 
described in [^], the linearized equations of motion describing gravitons in a general back- 
ground for quasi-topological gravity are fourth order in derivatives. As a result, one finds that 
the equations for the tensor perturbations in the black hole background now involve higher 
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powers of momentum, i.e., terms proportional to w^g^ and g^, as well as q'^d^. These quartic 
momentum terms will then dominate the large q limit. Indeed, in this case the full effective 
dispersion relation becomes 







2\f{p) + 2p\f'{p)) 



■ g2 (1 _ 2Xf{p) + 2\pf'{p) - VA/"(p)) 



-3puj 



f{p) ( f{p) - 2pf{p) + 27p'f"(p) + 48pV(')(p) + 12p^f^'\p) 



o4f(4)/ 



+3pV'(p) ( f'{p) + Qpf"{p) + 2p'f^'\p) ) + 6p'f"{p) 



+3/i 



M 

foo 



f{p) (f{p) - 2pf'{p) - 23pV"(p) - ^8p'f^'\p) - 12pV(')(p) 



+p'f'{p) ifip) - 2Apf"{p) - 12p'f^'\p) 



Joo 



-I2pp^i^q^{j'{p) + 2pnp)^ (c.^ 



(5.15) 



Hence as commented above, in the limit of large g, the contributions from the higher momen- 
tum terms, appearing in the last line above, come to dominate. In fact then, the dispersion 
relation reduces to that found for Einstein gravity ( ^.121) . Therefore we conclude that the 
dual excitations always propagate at less than the speed of light. A similar result was also 
found in |jl2| when considering a Weyl-tensor squared interaction added to the usual Einstein 
action. 

However, one might find this analysis somewhat suspect. In particular, this effective 
dispersion relation may not be well-defined here because of the very presence of these higher 
derivative terms which produced the simplification in the final step. Further in these higher 
derivative contributions, the asymptotic behaviour of the factor f"{r) + ^f'{r) ~ (ro/r)^ gives 
a very rapid decay and so perhaps modes propagating very close to the AdS boundary can 
evade our previous conclusion. 

To eliminate these potential loopholes, we now proceed with a more careful analysis by 
rewriting the equation of motion (5.11) in a Schrodinger form, following 15|. The first 
step is to isolate the contributions by rewriting eq. ( |5.11| ) as 



A{p, q2) 9Xp) + B{p, q2) 9p<A(/9) + C{p, q^) </>(p) + D{p, q^) w^^{p) = . 

Here we have defined the dimensionless frequency and momentum, 

U) q 
2^' ^~ 2^' 
Performing a change of coordinates and rescaling (pip) = Z{p)'ip(p) according to 

' 5pZ(p,q2) 



dy 
dp 



'D{py 



5pA(p,q2) 9pD(p,q2) B{p,q^) 



I A(p,q2) ' Z(p,q2) 
the eq. ( 5.16| ) takes the desired form 



4^(p,q2) 4Z)(p,q2) 2A(p,q2) 



(5.16) 
(5.17) 

(5.18) 



(5.19) 
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where o? = rx>^/q^. Note that m terms of the Schrodinger coordinate y, the horizon now 
appears at y — )• +00 and the asymptotic AdS boundary, at y = 0. In terms of the radial 
coordinate /), the effective potential can be expressed as 

"''^^'''^ Z)(p,q2) Z(p,q2) Z)(p,q2) Z(;,,q2) Dip,c^^) ' 
More concretely the effective potential is of the form, 

q^^(p,q^) = S|!!iM(5!)! (5.21) 

for some complicated functions rii and di. If we now take the large momentum limit we obtain 
the result 

= 4^ + 0(1/0- (5.22) 

loo 

Now this is precisely the effective potential which one would obtain for Einstein gravity in 
the large momentum limit. There is an infinite series of subleading corrections in 0(l/q^) 
which differ from Einstein theory but these terms are irrelevant in the limit of large q. This 
confirms our conclusion from the original analysis of the effective dispersion relation. 

Notice that in principle one has to worry that the boundary and large momentum limits 
do not commute. This is clear from the form of expression ( 5.2l| ), where the momentum 



appears in ratios. We will briefly comment on this in the discussion section. So perhaps a 
more subtle analysis may still find new constraints from demanding causality is respected in 
the dual CFT. 



Our discussion here has focused on the tensor modes ( 5.10 ). However, the subtleties 
regarding the boundary and large q limits carry over to the vector and scalar channels. In the 
large q limit, one obtains the same results as in the tensor channel. That is, in this limit the 
higher momentum terms dominate and one obtains a trivial dispersion relation, leading to 
no causality violation. Now one may wish to apply a more careful analysis for these modes, 
along the lines of that given above for the tensor channel. However, unfortunately for the 
vector and scalar channels, the previous analysis with an effective Schrodinger equation can 
not be applied in a straightforward way because of the appearance of higher powers of w. 

6. Holographic hydrodynamics 

In this section, we compute the ratio of the shear viscosity to entropy density for five- 
dimensional quasi-top ological gravity. By now, the holographic calculation of the shear viscos- 
ity is well understood. The first computations of this transport coefficient from an AdS/CFT 
perspective appeared in Q for Einstein gravity. These calculations were soon after extended 
to include higher curvature corrections to Einstein gravity, the first of example being the com- 
putation of the leading corrections to rj/s for the strongly coupled = 4 super- Yang- Mills 
theory |^. These computations were carried out for GB gravity [14, 17, 18 1 and also 
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higher order Lovelock theories |^T|, where the higher derivative terms need not be treated 
as smah corrections. Further investigations also provided increasingly efficient techniques for 
these calculations |]4l| , ^] In the following, we will use the 'pole method' of [^ ]. 

We begin with the metric for the planar AdS black hole given in eq. and which we 
write out again here 

ds^ = 4 dt^ + dx^ + dxl + dxl) + dr^ . (6.1) 

Recall that /(r) is determined by roots of the cubic equation in eq. (p.7|). Now it is convenient 
to transform to a radial coordinate z = 1 — r'^/r'^, with which the horizon is positioned at 
z = and the asymptotic boundary, at z = 1. The metric then becomes 

An important feature of these coordinates is that f{z) has a simple zero at the horizon. A 
Taylor expansion around z = yields 



where, e.g., /g = dzf\z=o- We present eq. ( |6.3D to establish a useful notation for the following. 
Following we perturb the metric (|6.2|) by shifting 



dx^dx + e e"*'^* dy , (6.4) 

where e is treated as an infinitesimal parameter. Then we evaluate the Lagrangian density, 
i.e., the entire integrand in eq. (|2.l| ) including ^/—g, on the shifted background to quadratic 



order in e. The presence of this off-shell perturbation (6^) produces a pole at z = in the 
(otherwise) on-shell action. The shear viscosity is then given by the 'time' formula [42| 

ri = —SttT lim ^^^^-0 ^ /-g ^-^ 

where Res2=o H denotes the residue of the pole in the Lagrangian density. Recall the Hawking 



temperature for the above black hole metric ( |6.2| ) is given in eq. (2.8). The final result of this 
calculation for quasi-topological gravity is 



?7 



1 - 2A/^ - 9^ [fi + 2/^'2 + 2/^ (/^' - 3/^)) J . (6.6) 

Now with the z coordinate, the cubic equation determining / can be written 

/(z)-A/(z)2-;./(z)3 = z(2-z) (6.7) 
and substituting in the Taylor expansion (|6.3| ), we can explicitly determine the coefficients: 
/^ = 2, /^ = -2(l-4A), = -24(A - - 2^) . (6.8) 
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The above expression {pM for the shear viscosity then becomes 



2B 



[l - 4A - 36/u(9 - 64A + 128X^ + 48^)] . 



(6.9) 



We readily verify that with ^ = 0, this result ( |6^ ) reduces to the expected result for GB 
gravity Q. 




Figure 2: Contours of constant rj/s shown in the allowed region of the gravitational couplings — see 
figure ^. The ratio increases in going from the lower-left to the upper-right in the plot. 



(6.10) 



Combining this result with eq. (|2.9| ), we find 



s 



^ [l - 4A - 36/i(9 - 64A + 128X^ + 48//)] . 



A contour plot of the ratio of the shear viscosity to the entropy density in space of gravitational 
couplings, A and /x, is shown in figure ^. From this plot, it is evident that rj/s is maximized 
in the upper-right corner of the the allowed region of couplings. This point corresponds to 



the intersection of the boundaries defined by the tensor and vector constraints, i.e., eqs. (5.5) 



and (5.6), respectively. At this point, one finds 

2838 , 246671 



2543 



A 



2684748 



6466849 
5714486118 



(6.11) 



Hence we find the minimum value for rj/s for the class of four-dimensional CFT's dual to 
quasi-top ological gravity is 



347182615788747017 1 , , , , 1 

~ (.4140) — . 

83858051009478068 1 47r ^ Mtt 



(6.12) 
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Note that this point is well away from the region of instabilities which will be discussed in 
the next section. Hence we expect that our calculation of r//s at this point is reliable. 

Within this class of CFT's, there is also a maximum value for r]/ s which appears to occur 



near the midpoint of the boundary produced by the scalar constraint (5.7). However, the 
point of the precise minimum lies in a region where, in the next section, we find that the 
uniform plasma is unstable — see figure ^. Hence our hydrodynamic calculations are not 
reliable at this precise point. Excluding the unstable region, it appears the maximum occurs 
very close to the point in GB gravity where rj/ s is maximized |17, That is, the maximum 
is near (A,/i) = (—1/8,0) where we find 



7/ 

s 



3 1 
2 4^ 



(6.13) 



7. Plasma Instabilities 

Even when the various consistency conditions of section ^ are satisfied, there remains the 
possibility that the black hole solution is unstable. The dual statement would be that an 
infinite uniform plasma is an unstable configuration for the CFT. Such an instability need 
not represent a fundamental pathology of the CFT but rather indicate that some interest- 
ing new dynamics arises in CFT plasma for certain values of the couplings. However, it is 
still important to identify such instabilities as they would invalidate the assumption of lo- 
cal thermodynamic equilibrium and for example, discredit the results of our hydrodynamic 
calculations in section |6|. 

The appearance of such instabilities for five-dimensional GB theory were first noted in 
[0, Although causality or positive flux constraints allow the GB coupling to be in the 
range —7/36 < A < 9/100, one finds that for A < —1/8 a new instability arises. Evidence 
for the latter was given as follows: First one writes the equation of motion for the tensor 
modes in an effective Schrodinger form, as was done in section above. For A < —1/8, the 
potential develops a small well where [/ < just in front of the horizon [i.e., near p = 1). For 
sufficiently large q, this well will support negative energy bound states which then correspond 



to unstable quasinormal modes, as described in [43|. Given that q is finite (and large), 
this instability indicates that the uniform plasma becomes unstable with respect to certain 
non-uniform perturbations. On the gravitational side then, this instability seems similar in 
certain respects to the Gregory-Lafiamme instability for black strings However, while 

the latter involves long wavelength modes, here the 'plasma instatiblity' occurs for arbitrarily 
short wavelengths. Examining the sound and shear channels, one finds that no additional 
instabilities arise in the consistent range, —7/36 < A < 9/100 |T^. The same analysis has also 
been extended to GB gravity in higher dimensions |l^. For D = 6, one finds similar range 
where the theory passes all the known consistency tests but the uniform plasma is unstable. 
However, for D > 7, all of the potential instabilities are pushed outside of the allowed range 
of the GB coupling. 
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In this section, we will provide a preliminary investigation of potential plasma instabilities 
for five-dimensional quasi-top ological gravity. Following the discussion above, our strategy 
will be to examine the tensor modes in detail using the effective Schrodinger equation, which 
was presented in eq. ( 5.1S| ). As described previously, because of the appearance of higher pow- 
ers of vo in the sound and shear mode equations, one cannot construct an effective Schrodinger 
problem in these cases. Hence a more elaborate analysis of the quasi-normal modes would be 
required to detect instabilities in these channels. 

We will separate our analysis into several different regimes, as we will find the behaviour 
of the theory will be quite different depending on the the sign of fi and the magnitude of 
q^. To see this, consider eq. ( 5.16 ). In particular, let us examine the coefficients ^(p, q^) and 
Z)(p,q2): 



dp 



1 - 2A/(p) - 3^/(p)2) 



d 



I2p^ {f{p) - 2pf'ip)) 

Joo 



A{p,c^' 



9p- 



P 



9/'(p) + 16p/(3)(p) + V/')(/o) 



(7.1) 



(7.2) 



The presence of the term proportional to fi q^ in the first expression creates the possibility that 
A (or D) may vanish, which we will see leads to a singularity in the Schrodinger potential. To 
see how this zero comes about, we first evaluate these functions at the horizon (i.e., p = 1). 
Using eq. ( ^.7)) , the polynomial defining f{p), and /? = Tq/f^, we find 



f{p):^-2ip-l)-il-4X)ip-lf + --- . 

With this result, we find the corresponding expansions of A and D at the horizon 

^0 



A~ Ao + O(p- 1), 



D 



4(p-l)2 



+ 0{p-l) 



(7.3) 



(7.4) 



with 



^0 = 24(3 - 8X)p ^ + [l - 4A - 36;u(9 - 64A + 128X'^ + 48p)] . 

Joo 

Hence we see that A vanishes at the horizon for a specific critical value of q^: 

2 _ 1 - 4A - 36^(9 - 64A + 128A^ + i8p) 
~ 24(3-8A)(-/i) ■ 

Now comparing with eq. ( |6.1C1| ), we see the expression in the numerator above is precisely 



(7.5) 



(7.6) 



Anrj/s. In the physically allowed region found in section 5.2, this ratio is always positive, as 
is the factor 3 — 8A — see figure |l]. Therefore we only have a valid solution for q^ when p is 
negative. 

Note that Aq vanishes when q = q^ and so from eq. ( [7. 4] ), both A and D vanish on the 
horizon at this point. For larger values of q, both A and D have a zero outside of the horizon 
{i.e., p < I) but the two zeros appear at different radii. 
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From the above discussion, it is clear that we should separate the analysis into three 
distinct regimes: i) > 0, ii) < and |q| < aiid iii) /x < and |q| > q^. 

7.1 Positive /i 

In the case where n > 0, the functions A and D are positive everywhere outside of the 
horizon. Hence, the effective Schrodinger potential ( 5.20| ) is well behaved everywhere in 



the range of interest, < /? < 1. To identify instabilities, it may seem that we can apply 
directly the strategy described above for GB theory of looking for a small negative dip in the 
effective potential just outside the horizon [|l4[ 15, |^]. However, there is a small subtlety 



which requires a more detailed investigation for quasi-topological gravity. For GB gravity, 
one further considers the limit of large q, which corresponds to the limit of /i — )• from 
the point of view of the effective Schrodinger problem. With this limit, any small negative 
dip will always lead to negative energy bound states as solutions to eq. ( ^.191 ) and hence 



the appearance of an instability [43|. However, as found in section [5.3| , the structure of the 
potential strongly depend on the value of the momentum for quasi-topological gravity. Indeed, 
our analysis showed that for sufficiently large momentum the effective potential reduces to 
that of Einstein gravity, i.e., f7 ~ f{p)/foo, and so no unstable modes would appear in this 
limit. That is, if one starts near an unstable point in GB theory but now with a small positive 
H, then generically the modes become more stable as the momentum is increased. Hence in 
general one will have to investigate the potential for finite values of q to find any unstable 
modes. 

Following the reasoning of we will still identify the unstable quasinormal modes as 
negative energy bound states in the effective Schrodinger problem. For a moment, let us 
consider applying the Bohr-Sommerfeld quantization rule for a zero-energy bound state: 

yjL./mi^)=(n-l).. (7.7) 

Here, we are assuming that the potential dips to below zero between the horizon at p = 1 and 
some lower turning point where U{p = po,'^'^) = 0. For the zero-energy state, the quantum 
number n would be a specific positive integer. More generally, if we were to evaluate the 
integral on the left-hand side, then the integer part of n on the right would count the number 
of negative energy bound states supported by the potential well. Hence our strategy here is 
to scan of the parameter space (A, n) for positive p. At each point, we vary q looking for a 
negative dip in the potential near horizon. If the dip becomes sufficiently deep to support a 



bound state {i.e., to satisfy eq. (7.7) with n = 1), we will take this to signal of an instability 
in the system. 

As discussed above for GB gravity, it was found that instabilities appear for A < —1/8. 
Hence we focussed our scan of parameters on small positive values of p in this regime and 
our numerical results are presented in figure ^. We should explicitly say that we expect that 
our Bohr-Sommerfeld analysis gives a good guide as to the unstable parameter space but it 
is difficult to assess how accurate the boundary of this region is in figure |3[ In this regard, it 
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is reassuring that the numerical curve reaches the A axis very close to A = —0.125 which, as 
indicated above, is where previous analysis [14, indicated that GB gravity should become 
unstable. 





Figure 3: Instability boundary for the positive fi estimated using the Bohr-Sommerfeld analysis. The 
(blue) region below the (red) curve is unstable. 



Note that the unstable region is very narrow just above the A axis but the height of this 
region increases as A becomes more negative. Intuitively, this behaviour arises because for 
GB gravity at ;U = 0, an infinitesimal negative dip first appears in the potential at A = —1/8 
and then becomes larger as A is decreased further. An infinitesimal dip appears is sufficient 
to support negative energy states in GB gravity because q can be taken arbitrarily large 
without effecting the effective potential. As discussed above, with nonvanishing /x, increasing 
the momentum makes the potential more stable or, in other words, decreases the size of dip 
in U. Hence one must balance this effect with the increase in the pre- factor q in eq. ( |7.7D . 
Therefore it is easier to find an instability for positive when when the initial size of the dip 
is larger at the corresponding point on the A axis. 



7.2 Negative fi and |q| < 

From our introductory discussion, we expect that the structure of the potential may be 
radically different for negative /i and particularly for |q| > q^. Here we will examine the 
approach to the critical momentum, q — )• q^, and the same Bohr-Sommerfeld analysis as 
above will show that unstable modes occur over a large part of this parameter regime. 

We are again looking for a negative potential well in front of the horizon, now for ^ < 
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and |ci| < q^. So to begin, consider the near horizon expansion of the effective potential ( ^.201) : 



U{p,q')^^{p-l) + —^{p-ir + --- , (7.8) 
^0 (^o)^ 



where Uq and Ui are some constants and Aq is given in eq. (7.5). An important point to 
notice is that as the momentum increases towards the critical value, — ^ 9c' have ^ 1 
and the above expansion becomes ill-defined. One must really perform a separate expansion 
for the special case q^ = qc, which yields 



Uip,<,l)^^^l^ + Oip-l). (7.9) 



where P and Q are polynomials in the couplings A and p. Their details are not important 
but for completeness we present them here 

P = 2(1 - AXfX + 3(-l + 4A) (-41 + 32A (19 - 76A + 96X^)) p 

-36(1393 + 12A(-1905 + 8A(1469 + 176A(-23 + 24X))))p^ 

-5184(161 + 24A(-49 + 96X))p^ - 1990656/ , (7.10) 
Q = 3(-l + 4A)(9 + 8A(-5 + 8X))p - 288 (20 + 3A (-189 + 16A (97 - 310A + 336A2))) p"^ 
-10368(-1 + 4A)(-45 + WiX)p^ + 248832/ . (7.11) 

The essential result here is that while for q^ < qc, the effective potential went to zero on the 
horizon, precisely at the critical momentum, it is tending to a constant value at the horizon. 
Hence the limit q^ — q^ is actually discontinuous. This unusual behaviour is illustrated with 
an example in figure 

The example in this figure also illustrates if this limiting constant value is negative, i.e., 
U{p = l,q^ = q^) < 0, then the effective potential develops a negative well in front of the 
horizon as we approach the critical momentum. Hence there is the possibility of developing 
instabilities in this regime. The dip in potential is largest with q^ = q^ and so we focus on 
this case. Once again, we will use the Bohr-Sommerfeld rule ( [7.7D to test for negative energy 
bound states. As before, we also have to worry that — )• as q^ — )• q^ and so we must 
perform a separate expansion to evaluate dy/dp at this point 



dy io , , . 3pQ{X,p) 

where Q was given above in eq. (|7ll| ) and 

R = 3pQ + 108 p{l- 4A)(3 - 8A) (3 - 80A(1 - 4A) + 320/ . (7.13) 
Hence the Bohr-Sommerfeld integral ( [7.7| ) yields an expression of the form 
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Figure 4: The effective Schrodinger potential in the tensor channel for /i — —0.0003, A —0.1 and 
different values of q^. Here — 7.2725. The effective potential at the horizon is tending to the 
constant value P/{(\lQ) ^ -0.797 

This integral produces logarithmic divergence indicating that there are an infinite number of 
unstable modes in this limit. Of course, this result is only true at strictly q = q^. However, 
by taking the momentum arbitrarily close to the critical value, the resulting negative dip in 
the potential will always support a large number of unstable modes as well. 

To summarize, we must exclude as unstable the region with /i < where the effective 
potential tends to a negative constant on the horizon at the critical value q,,. That is, the 
region where ;U < and P/Q < 0. The boundary of this region is well approximated by the 
curve 

2 ^ 14120 ,2 27472807424 . , 

~ A A-^. (7.15) 

^ 123 206763 28153056843 ^ ^ 

Putting this constraint for negative /i together with that obtained for positive in the previous 

section, one obtains the unstable region shown in figure ||. Here we are showing the region of 

instability together with the physically allowed region defined by the constraints in eqs. ( ^.5| - 

7.3 Negative /x and |q| > q^. 

With ^ < 0, we are of course free to take the momentum beyond the critical value. As noted 
above, when q^ > qci both A and D have a zero at some finite radius outside of the horizon 
but the two zeros appear at different radii. Let us label the two zeros as pA and po- In the 
following, we consider the case where pA > Pd, but in general this depends in detail on the 
values of A and p. However, the following analysis is easily adapted to the opposite situation 
where pA < Pd and the conclusions will be unchanged. 
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Figure 5: Unstable region (purple) for the tensor channel combining the results of sections [7.1| and 
7.2, superimposed onto the region (pink) allowed by positivity of energy flux. In the section [77^ , the 
analysis is extended to |ci| > q^, and it appears instabilities are present throughout the entire /i < 
region. 



Let us first consider the effective potential in tlie vicinity of Z)'s zero. Expanding the 



various coefficients in eq. (5.16) around p ~ pd gives 



A = Aq + 0{p- pd) , 
B = Bo + O{p- pd) , 
C = Co + 0{p-pd) , 
D = -Di{p-pn) + 0{ip-pD?) . 



(7.16) 



It is important to note that the constants Aq and Di are both positive, so that dy/dp in 
eq. ( |5.18D is well defined in the region p < pr,. Using the expressions above, along with 
eq. (|5.18|), we find 



VD-y 
Z 




3 V ^0 

{pd-pT'I'. 



iPD - pf' 



(7.17) 
(7.18) 



Finally, with eq. ( 5.2(]| ), we obtain the effective potential 



^^U{p) 



5/16 



Di {p-poT 
5/36 

iVD - vY ' 



(7.19) 
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Hence we find that the effective potential contains a singularity at p = po but the structure of 
the singularity is surprisingly simple. In particular, the coefficient is completely independent 
of the parameters, A, p and q, but of course, the latter still determine the precise location of 
the singularity. 

We can perform a similar analysis for the effective potential in the vicinity of the zero in 
A. We begin by expanding around p pA 

A = -A,{p-pA) + 0{{p-pDf) , 
B = -Ai + 0{p- pa) , 
C = Cq + 0{p-pa) , 

D = -Do + O{p- pa) . (7.20) 
for some positive constants Dq and Ai. Note that the leading coefficient in B has the special 



form Bq = Ai, which follows from the original equation (5.11). Then in this case 



yA^2^^^ip-PA)'^\ (7.21) 
Zc^ip-PA)-'/' (7.22) 



and the potential becomes 



Do p- PA 
1/4 

(y - VA 



\2 



(7.23) 



Here again we find that the effective potential contains a singularity with a very simple form 
at the zero of A. The parameters, A, p and q, fix the precise location of the singularity but 
the overall coefficient is independent of these. 

Implicitly the above analysis assumes that p > pA where both A and D are negative. 
Similarly in consider the zero in D, we assumed p < po where both A and D are positive. 
Special consideration must be given to the range Pd < P < PA where D < and A > 0. 
To accommodate the latter, we must modify the construction of the effective Schrodinger 
equation ( |5.19| ) slightly. In particular, in eq. ( ^.18| ), the definition of the new coordinate is 
replaced by 

dy_ I D{p,^^) 

Tp-^-^iiW)- ^'-''^ 

The final result can then be most simply expressed as flipping the sign of both the effective 
potential and the effective energy of eq. ( ^.19| ). That is, in this region, the Schrodinger 
equation becomes 

--^5^'V'(2/)+ [-C/(2/,c,2)] V^(y) = [-a2] ^(y) , (7.25) 
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where U is given by precisely the same expression in eq. ( 5.20 ) and = t»^/q^, as before. 



For present purposes, the behaviour near the zeros p = and pA are of primary interest. 
Carefully keeping track of the signs, one finds that near these points, the Schrodinger equation 
( 7.25| ) has precisely the same form as above, with a singular and attractive potential. That 
is, in the vicinity of these zeros, we may write the Schrodinger equation as 



(7.26) 



where 7 = 5/36 for p = po and 1/4 for p = pA- In either case we have shifted the y coordinate 
to put the singularity at the origin. The 1/y^ potential has been studied extensively in the 
literature |45]. Remarkably for a one-dimensional Schrodinger equation, the value 7 = 1/4 
marks the boundary between the conformal regime when the potential is repulsive or weakly 
attractive and the regime with a sufficiently attractive potential where discrete bound states 
appear and conformality is lost. 

The key point in presenting eq. ( 7.26| ) that this equation applies here for both positive 
and negative y. However, in this one-dimensional setting, the wave-function must propagate 
through the singularity at y = 0. So one should approach the problem by solving for y < 
and y > independently and then matching the solutions with an appropriate boundary 
condition at the origin. This procedure is most easily demonstrated for the zero in D, in 
which case 7 = 5/36. Then eq. (7.26) has two independent solutions 



V'(y) 



di{-y)'/' + d2{-yf/' 



for y > , 
for y < . 



Using eqs. (7.17) and ( |7.18 ) above, we can convert these back to the radial profile of the 
original tensor mode, i.e., (j) = Zip, in the vicinity p ~ po 



4>{p) = ^1 + 52{p - Pd) 
= 5i + 52{p - pd) 



for p < Pd, 
for p> PD ■ 



(7.27) 



Matching of solutions at p^i is achieved straightforwardly by imposing continuity of the radial 
profile and its first derivative, i.e., 5i = 5i and 82 = 52- Hence while some care must be taken, 
the zero in D presents no real difficulty in finding well-behaved solutions. 

Next we turn to the zero in A, where the situation is more subtle. In this case, 7 = 1/4 
in eq. ( [7.26| ) and the two independent solutions are 

ai y^'^ + a2 log(y) 
ai 



V'(y) 

Hv) = ai i-y)^'^ + 02 i-yf'^ log(-y) 



for y > , 
for y < . 



Using eqs. (|7!2l| ) and (^), these expressions are translated to the radial profile 
the vicinity p ~ 



<l){p) = ai + 02 log(p - pa), 
= Oil + a2log(pA - p), 



for p> PA, 
for p < PA- 



Z-p in 



(7. 
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The matching at pA is shghtly more involved because of the logarithmic behaviour of these 
solutions. Integrating the equation of motion ( ^.16| ) around pA, we obtain 

- PA)dp^\'% = dp (Doto^-Co)cl). (7.29) 

As e — )• 0, the right-hand side vanishes and the vanishing of the left-hand side requires 
a2 + 02 = 0. Another natural boundary condition is that the flux of probability in the 
effective Schrodinger problem should be continuous as the wave-function propagates through 
the singularity at ?/ = 0. If we express this flux in terms of the radial profile and the 
coordinate p, then we require 



lm[<p*A,{p-pA)d,cf\2-^_l = 0. 



Continuity then leads to Im(a;Ja2) = Im(Q!jQ;2) and so it seems natural to set ai = — di as 
well.^ Our final solution in the vicinity of p ~ then takes the form 

(l)= [ai + 02 log \p - pa\) sgn(p - pa) ■ (7.30) 

Again the constants ai and 02 are arbitrary but our analysis shows appropriate boundary 
conditions at this point will produce a suitable physical solution. 

At this stage, we have shown that despite of the appearance of new singularities in the 



equation of motion (5.16) in this regime {i.e., p < and |q| > q^), we may more or less 
straightforwardly solve for the radial profile (p. The precise solutions and the corresponding 
quasinormal eigenfrequencies to will still be set by the boundary conditions on the field at 
the black hole horizon p = 1 and at the asymptotic boundary p = 0. To better understand 
these boundary conditions, we now return to the overall behaviour of the effective Schrodinger 
potential. 

Figure ^ provides an example of the effective potential in the desired regime. One point 
which the figure illustrates is that for large momenta ^ q^, the general discussion of 



section |5.3| still applies here and over most of the range the effective potential approaches 
that of Einstein gravity, i.e., U ~ f{p)/ foo- However, the figure also makes evident the 
singularities extensively discussed above, which appear as sharp dip ai p = pA and po- As is 
typical of q^ ^ 9ci zeros are very close together and in fact it is difficult to resolve the 
two distinct singularities in the example given in figure |6|. Intuitively, one expects that this 
deviation of the Einstein potential will provide a small perturbation and so there will be a set 
of stable modes whose quasinormal frequencies and radial profiles deviate only slightly from 
the solutions in Einstein gravity. 

However, the singularities introduce a new boundary surface into the problem and we 
argue that this also leads to a additional set of new unstable modes, as follows: The Einstein 



^This is not a unique solution for this constraint and so it may be that a further boundary condition should 
be applied to single out this result as the unique physical solution. 
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Figure 6: Comparison of the effective potential for /i < and > (red) with that for Einstein 
gravity (black). The red curve is given for A = —0.1, /i = —0.0003 and q = 40. For these gravitational 
couplings, the critical momentum is approximately q,, ~ 7.27. The inset shows a close-up of the 
effective potential to resolve the two separate singularities at pA = 0.328378 and po = 0.327955. 



(Note that in keeping with the discussion at eq. ( |7.25D , the sign of U in the inset has been reversed 
on the interval pA > p > Pd-) 



potential vanishes at the horizon and so solving the Schrodinger equation ( 5.19| ) with a neg- 
ative energy yields two independent solutions, one which grows exponentially (and diverges 
at the horizon) and another which decays. Similarly, there are two asymptotic solutions, 
one which grows and another which decays as one approaches the AdS boundary. Now the 
potential is smooth throughout < p < 1 if we are considering pure Einstein gravity. As a 
result, the solution which decays at the horizon is precisely that which grows at the asymp- 
totic boundary and vice versa. Hence one finds no normalizable solutions with a negative 
effective energy which agrees with the result that the black hole is stable in Einstein gravity. 
However, if we consider quasi-topological gravity with p < and |q| > q^i while the effective 
potential is well approximated by the Einstein potential for most radii, a dramatic difference 
arises with the appearance of the singularities at p = pA and po- Now we have the possi- 
bility of matching a solution in the range p > pA which decays at the horizon to a solution 
in the range p < pD which decays at the asymptotic boundary. Because of the sign of the 
effective energy in the Schrodinger equation ( 7.25 ) is flipped in the range pD < p < PA, the 
solution would be oscillatory in this interval and tuning should allow us to match onto the 
decaying solutions with the boundary conditions indicated above. Hence we should also be 
able to construct an infinite set of negative energy states which are localized near the interval 
Pd < P < PA- It appears that these unstable modes would naturally be regarded as the 
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progeny of the diverging number of unstable modes which were found to accumulate in the 
limit — ^ 9c' ™ previous section. In any event, our conclusion is that generally when fj, 
is negative, the black holes in quasitopological gravity have unstable quasinormal modes with 
|q| > q^. 

To summarize, our analysis in this section indicates that instabilities appear in the tensor 
channel throughout the lower half of the space of couplings with /x < and in the narrow 
sliver shown in figure ^ with // > 0. However, we should re-iterate that our analysis here only 
represents a preliminary investigation of the potential plasma instabilities in five-dimensional 
quasi-topological gravity. A thorough and detailed analysis of the quasinormal modes is 
required to validate the results derived here and to develop a full picture of the unstable 
modes. 



8. Discussion 



In this paper, we have begun an examination of the holographic properties of quasi-topological 
gravity. The main new feature of this toy model is that it allows us to examine CFT's in 
which the flux parameter in eq. ( |4.lD is nonvanishing. In this case, the dual CFT cannot 



be super symmetric |11] and so these new models allow us to begin exploring holography in a 
context which is fundamentally nonsuper symmetric. In this regard, quasi-topological gravity 
differs from Lovelock theories which are consistent with supersymmetry [|^. Of course, 
conformal fixed points are believed to occur for a wide variety of nonsupersymmetric gauge 
theories [^] . Further one might speculate that extending some of these to large Nc and strong 
coupling may generate a holographic dual close to Einstein gravity [48|. Hence, in the spirit 
of exploring quasi-topological gravity as a toy model, one may gain new insights into such 
conformal fixed points. 

One aspect of physics which we explored was the hydrodynamic properties of the CFT 
plasma. In particular, we calculated the ratio of shear viscosity to entropy density ( |6.10D . 
Similar results have already been found in a framework where the curvature-cubed terms 
were treated perturbatively [Q. Eq. (|6.10| ) should reduce to these results in the regime of 
small couplings where one only keeps the terms linear in A and /i. Our full nonperturbative 
calculation is straightforward using the techniques developed in and produces a final 
expression which contains contributions which are nonlinear in the gravitational couplings, 
e.g., proportional to /iA^. This differs from cubic Lovelock theories in higher dimensions 
which also contain curvature cubed interactions. In fact, in this case, rj/s is independent of 
the coupling constant controlling the curvature-cubed interactions and remains simply 

linear in the curvature squared coupling, as in GB gravity p^ , 16, 1^. Of course, these theories 
are also distinguished from the present case since t4 remains zero in the Lovelock theories 
despite the appearance of the curvature cubed interaction [^] . 

Considering the value of rj/s when t4 ^ 0, we find no dramatic behaviour. The value 
smoothly increases or decreases as we move away from the axis into the (A, fj.) plane, as 
illustrated in figure |2[ One point that this analysis makes clear is that rj/s = 1/Att is simply 
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a codimension one surface (here a contour) which cuts through the space of gravitational 
couphngs or alternatively through the space of parameters which differentiate the dual CFT's. 
The main feature that distinguishes this surface is that it runs through the point where the 
bulk theory is Einstein gravity, which we favour as theorists. This illustrates that even if r//s 
was found to be precisely 1/4tt for some system arose in nature {e.g., the quark-gluon plasma 
or a trapped atomic gas), there is no guarantee that the holographic dual would anywhere be 
close to a theory of Einstein gravity. 

This holographic model also illustrates the point that the CFT plasma can readily achieve 
r]/s < 1/4:11. Of course, even though the originally conjectured KSS bound on ij/s has been 
proven incorrect, there are still general arguments to suggest that this ratio should satisfy 
some lower bound |3|, Hence the question naturally arises as to the precise nature of 
such a bound. Holographic models provide an excellent theoretical framework to study this 
question, as they provide access to a variety of strongly coupled fluids in this 'KSS regime' 
where rj/s is unusually small. 

In the present case of holographic fluids modeled by quasi-topological gravity, this ratio 
reaches a minimum value at the upper corner of the allowed parameter space with r//s ~ 
(0.4140)/(47r), as given in eq. ( |6.12 ). It may be of interest to note that at the point, the 
dual CFT has t2 = while ^4 reaches its maximal value with = 15/2. It is clear that 
this value does not represent a fundamental bound. If one explores holographic models 
with GB gravity in higher dimensions, one finds the minimum ratio is r]/s ~ (0.4139)/(47r) 
for D ~ 9.207 |16, 17|. This analysis was also extended to Lovelock gravity with a cubic 
curvature interaction for D > 7 |^l|. Initially here, it appeared that rj/s could be pushed to 
zero (or even negative values) but taking care to account for plasma instabilities, one finds 
that the hyrdodynamic results are only reliable down to r//s ~ (0.3938)/(47r) |5^. The cubic 
interaction of quasi-topological gravity can also be extended to higher D > 7 [|^] and so it 
would be interesting to examine the contributions of such interactions to higher dimensional 
holographic models. One might note the Lovelock models introduce additional parameters to 
distinguish the dual CFT's but still have = 0. While it is perhaps not surprising then, these 
results explicitly demonstrate that many more parameters in the CFT will effect the value of 
r]/s than simply A, B and C, the three couplings which fix the the three-point function of the 
stress tensor. 

While the various holographic models above all seem to point to a minimum value around 
r//s ~ (0.4)/(47r), it seems clear that this collection of models only explores a limited param- 
eter space. One might imagine that one can continue to systematically lower rj/s by contin- 
uing explore a wider range of CFT's by adding more and more interesting couplings in the 
dual gravitational theory (and without introducing any other pathologies in the holographic 
model). In fact, there seems to be some evidence in favour of such a scenario at least in very 
high dimensions | ]51| . 

There is one difference between quasi-topological gravity and GB gravity (or the more 
general Lovelock theories) which is particularly striking. All of the Lovelock gravity models 
are distinguished by having second order equations of motion while the general equations 
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in quasi-topological gravity are fourth order. It seems that the latter is inevitable in order 
to produce a holographic theory where t4 7^ since the Lovelock theories are in fact the 
most general gravitational theories with second order equations of motion fS^ . To develop 
some intuition for such higher order equations, we might establish an analogy with a higher- 
derivative scalar field equation (in flat space) 

Here we imagine is some high energy scale and a is the dimensionless coupling that controls 
the strength of the higher-derivative term. The (flat space) propagator for this scalar can be 
written as 

111 , , 

(8.2) 



g2(l-ag2/M2) g2 q2_M^/a' 

Now the l/g2 pole is associated with the regular modes which are easily excited at low 
energies. The second pole l/(g^ — M^/a) is associated with ghost modes that appear out at 
the high energy scale. Depending on the sign of a, these new modes may have a regular mass 
(a < 0) or be tachyonic (a > 0). Further writing out the dispersion relation for the ghost 
modes, we have 

g2-M2/o = -w2^(P)2-MVa = 0. (8.3) 

As noted, when a is negative, the mass above has the 'right sign' and these modes only go 
on-shell when cj^ ~ M'^/\a\, i.e., at very high energies. On the other hand, if a is positive, 
the modes are 'unstable' and in this case, we can bring these modes on-shell above a certain 
threshold of large spatial momentum, i.e., (fc*)^ ~ M'^/a. Comparing this discussion to our 
analysis of the tensor channel equation for quasi-topological gravity, the coupling // for the 
curvature-cubed interactions would play a role analogous to a above. In parallel with the 
present scalar theory in section |7.3| , we found that a new set of unstable modes appears above 
a certain momentum threshold for a particular sign of /i, i.e., > for fi < 0. 

One important point that arises in the scalar field model is that the extra high energy 
modes are ghosts for either sign of a, as seen from the overall sign of their contribution to 
the propagator ( |8.2D . Hence a natural worry would be that ghosts must also appear in quasi- 
topological gravity but we will argue that this is not the case, in the following. In section 0, 
we have found new unstable modes in quasi-topological theory. This certainly indicates that 
working on the uniform black hole background is problematic but it is not clear that they 
indicate that there is a fundamental pathology in the form of ghost modes. One important 
difference between the equations in quasi-topological gravity and in the simple scalar model 
are that the former are not Lorentz invariant (in the gauge theory directions). Of course, 
the lack of Lorentz invariance is not a surprise since the black hole background is dual to a 
uniform finite temperature plasma which defines a preferred reference frame. Hence the same 
statement would apply even if we were considering graviton modes in a black hole solution of 
Einstein gravity. However, the discrepancy is more significant for the higher derivative terms 
here. For example, while a term appears in eq. (|5l^ ) there is no m contribution. Hence 
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it is not clear whether the additional instabilities appearing in the black hole background are 
associated with ghost modes. 

Let us consider this point further. At zero temperature, the background spacetime re- 



duces to simply AdSs and as noted before with eq. ( 2.10 ), the linearized equations of motion 
reduce to the second order equations of Einstein gravity. Hence in this limit, both the higher 
derivative contributions and the additional unstable modes vanish for any values of A and 
H- As noted in |25[ |, the higher derivative terms appear in the linearized equations of motion 
through couplings to the background curvature when one is considering fluctuations around a 
nontrivial background spacetime. In particular, these terms arise from a nontrivial Weyl cur- 
vature in the background, e.g., for a transverse traceless mode {i.e., V^hab = 0) and h^a = 0), 



the four-derivative terms can be written as |25] 



.6 



C^'^ hde;cfiab) + 2 ''''q,de\b) + 2 O^h'^'^C^adb + gab (□/icd);e/ C"^''^ • (8-4) 

Hence the specific features of any instabilities will always depend on the details of the back- 
ground geometry under study. It may be of interest to explicitly repeat the analysis of section 
for other nontrivial configurations, i.e., one simple example would be the confining phase 



represented by the AdS soliton [53 



Another interesting aspect of these higher derivative terms is that they will vanish 
in the asymptotic region with AdS boundary conditions, since the Weyl curvature will vanish 
there. Hence any unstable modes associated with these terms will be 'confined' to the interior, 
e.g., near the horizon. Hence in the dual CFT, the instabilities will be associated with 
dynamics of infrared excitations and will be insensitive to the ultraviolet details of the theory. 
Given this observation and the previous discussion, it seems then that these problems cannot 
represent a fundamental pathology, i.e., ghosts, in the theory. Hence we would conclude that 
the analogy with the scalar field in eq. (8.1) is simply deficient in this respect. 



Still the higher derivative terms and the resulting instabilities are a worrying aspect of 
quasi-topological gravity. In particular, we have made a preliminary analysis of the shear and 
sound channel equations of motion. In this case, it appears that the coefficient analogous A 
in eq. ( |5.16| ) can again have go to zero but when fj, is positive. Hence one would be tempted 
to conclude that instabilities will now appear for positive fi. Unfortunately, this would mean 
that the plasma is only stable with fj, = 0. This makes the need for a detailed analysis of the 
quasinormal modes even more pressing to develop a full picture of the instabilities. 

To close, let us re-iterate that quasi-topological gravity was introduced as a toy model to 
study extensions of the usual AdS/CFT correspondence. We have not identified an approach 
by which the new gravitational action ( |2.1D emerges from a UV complete theory and we 
have no reason to expect that the new theory is radiatively stable. If one were tempted to 
construct a full quantum version of quasi-topological gravity, another issue which would need 
to be addressed is the appearance of several different vacua (corresponding to the roots of 
eq. (^.5|)) and in particular vacua in which the graviton is a ghost. It is not clear what the 



^Here we adopt the standard notation: T^ab) = | (Tat + Tba). 
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role of these vacua would be in, e.g., a path integral formulation of the theory. All of these 
considerations as well as the appearance of instabilities at large momenta reinforce the idea 
that quasi-topological gravity (as well as Lovelock theories of gravity) should only be treated 
as toy models which may give us insight into the long wavelength physics, e.g., hydrodynamic 
behaviour, of strongly coupled conformal fixed points. 
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